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Abstract 

We present the DTFE public software, a code for reconstructing fields from a discrete 
set of samples/measurements using the maximum of information contained in the point 
distribution. The code is written in C++ using the CG Alj^ library and is parallelized 
using OpenMP. The software was designed for the analysis of cosmological data but can 
be used in other fields where one must interpolate quantities given at a discrete point 
set. The software comes with a wide suite of options to facilitate the analysis of 2- 
and 3-dimensional data and of both numerical simulations and galaxy redshift surveys. 
For comparison purposes, the code also implements the TSC and SPH grid interpolation 
methods. The code comes with an extensive user guide detailing the program options, 
examples and the inner workings of the code. The DTFE public software and further 
information can be found at |http : //www . astro . rug . iil/~voroiioi/DTFE/dtf e . htmll 
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1 Introduction: the DTFE method 



The Delaunay Tessellation Field Estimator (from now on DTFE) (see Schaap & van de Wey- 
gaert|2000 van de Weygaert fc Schaap|2009 Cautun et al.|2011 ) represents the natural method 
of reconstructing from a discrete set of samples/measurements a volume-covering and contin- 
uous density and intensity fields using the maximum of information contained in the point 
distribution. This interpolation method works for any scalar or vector fields (e.g. velocity, 
temperature, pressure) that are defined at the positions of a discrete point set. Moreover, 
if the point distribution traces the underlying density field, DTFE offers a local method in 
determining the density at each point and accordingly also in the whole volume. 
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^Computational Geometry Algorithms Library, http://www.cgal.org 
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Figure 1: An illustration of the 2D Delaunay tessellation of a set of particles from a cosmological 



numerical simulation. Courtesy: Schaap (2007). 



The DTFE method was first developed by Willem Schaap and Rien van de Weygaert (Schaap 



&: van de Weygaert||2000 ) to be used on various astrophysical applications, but can also be used 
in other fields where one must interpolate using quantities given at a discrete point set. The 
DTFE method is especially suitable for astrophysical data due to the following reasons: 

1. Preserves the multi-scale character of the point distribution. This is the case in numerical 
simulations of large scale structure (from now on LSS) where the density varies over more 
than 6 orders of magnitude and for a lesser extent for galaxy redshift surveys. 

2. Preserves the local geometry of the point distribution. This is important in recovering 
sharp and anisotropic features like the different components of the LSS (i.e. clusters, 
filaments, walls and voids). 

3. DTFE does not depend on user defined parameters or choices. 

4. The interpolated fields are volume weighted (versus mass weighted quantities in most 
other interpolation schemes). This can have a significant effect especially when comparing 



with analytical predictions which are volume weighted (see Bernardeau & van de Weygaert 
T996|. 
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The first two points can be easily seen in Fig. [T| where starting from the point distribution 
(upper left insert), one constructs the Delaunay tessellation (upper right insert) and than zooms 
on a high density region (lower right and left inserts). 



DTFE density SPH density 




Figure 2: A comparison between the DTFE (left panel) and SPH (right panel) densities through 
a 80 X 80 X 2 Mpc/h slice in an N-body simulation. Both density fields are represented on a 
logarithmic scale. The SPH density was computed using an adaptive smoothing length that 
involves the closest 40 neighbors. 



2 The DTFE public software 

The code accompanying this document is a C++ implementation of the DTFE method of in- 
terpolating from a discrete point set (in 2 or 3 dimensions) to a grid. The software was designed 
using a modular philosophy and with a wide set of features that can easily be selected using the 
different program options. These features were constructed to analyze the data using multiple 
methods and to manipulate and split the data such that it can be useful in dealing with a wide 
variety of problems. The DTFE code is also written using OpenMP directives which allow it 
to run in parallel on shared-memory architectures. The software can be used as a standalone 
program or as an external library. 

The code comes with complete documentation and with a multitude of examples that detail 
the program features. Moreover, a help desk is available for information and assistance for 
troubleshooting problems. The DTFE software can be downloaded from http : // www . astro .J 
[rug ■ nl /~voronoi /DTFE/ dtf e . html, The same address contains the user guide as well as help 
and contact forms. 
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Figure 3: A map of the DTFE computed velocity flow (left panel) and velocity divergence (right 
panel) corresponding to the density field shown in Fig. [2] The left panel shows the velocity 
vectors colored according to the velocity magnitude. The right panel show the absolute value 
of the velocity divergence such that it can be shown on a logarithmic scale. 



3 DTFE software features 

The DTFE software has a wide range of features designed to facilitate the work with cosmo- 
logical data, both from numerical simulations and from observations: 

• Works in both 2 and 3 spatial dimensions. 

• Interpolates the fields to three different types of grids: 

1. Regular rectangular and cuboid grid — useful for numerical simulation. 

2. Redshift cone (spherical coordinates) grid — useful for galaxy redshift survey or for 
simulating observations. 

3. User given sampling points — can describe any complex or non- regular sampling 
geometry 

• Uses the point distribution itself to compute the density and interpolates the result to 
grid. 

• Each sample point has a weight associated to it to represent multiple resolution N-body 
simulations and observational biases for galaxy redshift surveys. 

• Interpolates the velocity, velocity gradient, velocity divergence, velocity shear and velocity 
vorticity (see Fig. Islfor an example of the DTFE velocity and velocity divergence fields). 
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• Interpolates any additional number of fields and their gradients to grid. 



• Periodic boundary conditions. 

• Zoom in option for regions of interest. 



• Splitting the full data in smaller computational chunks when dealing with limited CPU 
resources. 



• The computation can be distributed in parallel on shared-memory architectures. 



For comparison purposes, the software comes also with the TSC(Hockney & Eastwood 
1981) and the SPH(Monaghan 1992) interpolation methods (see Fig. |2jfor a comparison 



of DTFE and SPH density fields). 



• Can return the Delaunay tessellation of the given point set. 

• Easy change of input /output data format. 

• Easy to use as an external library. 

• Extensive documentation of each feature. 



4 DTFE software license 



DTFE is free software, distributed under the GNU General Public License This implies that 



you may freely distribute and copy the software. You may also modify it as you wish, and 
distribute these modified versions as long as you indicate prominently any changes you made 
in the original code, and as long as you leave the copyright notices, and the no-warranty notice 
intact. Please read the General Public License for more details. Note that the authors retain 
their copyright on the code. 



If you use the DTFE software for scientific work, we kindly ask you to reference the DTFE 
method and code papers: Schaap & van de Weygaert (2000), van de Weygaert & Schaap (2009) 
and Cautun et al. (2011). 
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